A novel method to create a vortex in a Bose-Einstein condensate 
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It has been shown that a vortex in a BEC with spin degrees of freedom can be created by 
manipulating with external magnetic fields. In the previous work, an optical plug along the vortex 
axis has been introduced to avoid Majorana flips, which take place when the external magnetic field 
vanishes along the vortex axis while it is created. In the present work, in contrast, we study the 
same scenario without introducing the optical plug. The magnetic field vanishes only in the center 
of the vortex at a certain moment of the evolution and hence we expect that the system will lose 
only a fraction of the atoms by Majorana flips even in the absence of an optical plug. Our conjecture 
is justified by numerically solving the Gross-Pitaevskii equation, where the full spinor degrees of 
freedom of the order parameter are properly taken into account. A significant simplification of the 
experimental realization of the scenario is attained by the omission of the optical plug. 
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I. INTRODUCTION 

Alkali atoms become a superfiuid upon Bose-Einstein 
condensation (BEC) [1, 2]. The superfiuid properties of 
the system were considered to be essentially the same 
as those of superfiuid 4 He, in spite of the fact that the 
former is a weakly-coupled system while the latter is cou- 
pled strongly. In contrast with 4 He, however, alkali atoms 
have internal degrees of freedom attributed to the hyper- 
fine spin F and, accordingly, the order parameter has 
2\F\ + 1 components [3, 4]. The atoms 23 Na and 87 Rb 
have = 1, for example, and their order parameters 
have three components, similarly to the orbital or the 
spin part of the superfiuid 3 He. 

These degrees of freedom bring about remarkable dif- 
ferences between the BEC of alkali atoms and that of 4 He. 
The hyperfine spin freezes along the direction of the lo- 
cal magnetic field when a BEC is magnetically trapped. 
If a BEC is trapped optically, on the other hand, these 
degrees of freedom manifest themselves and various new 
phenomena, such as the phase separation between the 
different spin states that have never been seen in super- 
fluid 4 He, may be observed. 

It was suggested in [5] and [6] that a totally new pro- 
cess for formation of a vortex in a BEC of alkali atoms is 
possible by making use of the hyperfine degrees of free- 
dom to "control" the BEC. Suppose a BEC is confined 
in a Ioffe-Pritchard trap. Then a vortex state with two 
units of circulation can be continuously created from the 
vortex free superfiuid state by simply reversing the axial 
magnetic field B z , parallel to the Ioffe bars. This field 
is created by a set of pinch coils and can be easily con- 
trolled. Four groups have already reported the formation 
of vortices with three independent methods [7-11]. Our 
method is totally different from the previous ones in that 



the hyperfine spin degrees of freedom have been fully uti- 
lized. The present paper is a sequel to [6]. Here, we con- 
duct further investigations on the formation of a vortex 
by reversing the axial field B z . An optical plug along the 
axis was introduced in [5] and [6] to avoid possible Ma- 
jorana flips near the axis, that may take place when B z 
passes through zero. It is expected, however, that trap 
loss due to Majorana flips may not be very significant if 
the "dangerous point" \B\ = is passed fast enough and 
that considerable amount of the BEC remains in the trap 
as the result. This process should not be too fast, how- 
ever, such that the adiabatic condition is still satisfied. 
It is a very difficult task to introduce a sharply forcused 
optical plug along the center of the condensate whose ra- 
dius is of the order of a few microns. Accordingly we 
expect that experimental realization of our scenario will 
be much easier without the optical plug. In the present 
paper, therefore, we analyze our scenario by numerically 
integrating the multi-component Gross-Pitaevskii equa- 
tion. 



This paper is organized as follows. In Sec. II, we out- 
line the order parameter and the Gross-Pitaevskii equa- 
tion for an |F| = 1 BEC. In Sec. Ill, the Gross-Pitaevskii 
equation is integrated numerically to analyze the time de- 
pendence of the order parameter. It will be shown that 
merely half of the condensate is lost from the trap if the 
time dependence of the external magnetic field is chosen 
properly. Sec. IV is devoted to conclusions and discus- 
sions. In the Appendix, it is shown that the formation 
of a vortex in the present scenario may be understood in 
terms of the Berry phase. 
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II. ORDER PARAMETER AND 
GROSS-PITAEVSKII EQUATION OF \F\ = 1 BEC 

Let us briefly summarize the order parameter and the 
Gross-Pitaevskii equation for a BEC with the hyperfme 
spin F = 1 to make this paper self-contained. The read- 
ers should be referred to [3] and [6] for further details. 

The order parameter of a BEC with |F| = 1 has three 
components ^ (i = — 1,0, +1) with respect to the basis 
vectors \i) defined by F z \i) — i\i). The order parameter 
|\&) is then expanded as 

l*>= E **!»)• 

i=0,±l 

It turns out, however, that another set of basis vectors 
| a) (a = x,y,z), defined by F a \a) = 0, is more conve- 
nient for certain purposes. The order parameter is now 
expressed as 



i*>= E *<» 



a— x,y ,z 

The transformation matrix from {^i} to {^ a } is found 
in [6]. 

The most general form of the Hamiltonian for \F\ = 1 
atoms, that is rotationally invariant except for the Zee- 
man term, is 



H - fiN 



1>l 



.^V 2 - M Ua + f(^) 2 
Zm I I 



92 



\lpai>a\ 2 + ie abc UJ Lc 4>llp b 



<Fr, 



(1) 



where hui La = j^Ba, 7 M ~ Ms/2 being the gyromagnetic 
ratio of the atom. The coupling constants are given by 



4irh 2 



9i = 



-a 2 , 92 



4irh 2 clq — ai 



(2) 



where a 2 = 2.75 nm and eto = 2.46 nm [12]. 

The Heisenberg equation of motion derived from the 
above Hamiltonian is 

\ 2m J 

+92('ipb^b)lpi+i£abc^Lci>b- (3) 

By taking the expectation value of the above equation 
in the mean-field approximation, (ipipip) ~ ^I'^E'^', where 
*a = {%), we obtain 



ih 



dt 



fir 2 (*) 2 *a + ie a bcULc^b- 



(4) 



This is the fundamental equation which we will use in 
the rest of this paper. 
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FIG. 1: The schematic time dependence of the axial field B z . 
The field B z reverses the sign in an interval of time T, while 
the quadrupole field B± = B'r remains fixed. The field B z (t) 
is fixed at the value B X (T) = -B z (0) for t > T for further 
evolution. 



III. CREATION OF A VORTEX 
A. Magnetic Fields 

Suppose we confine a BEC in a Ioffe-Pritchard trap 
with a quadrupole magnetic field 



B ± (r) = 



B±(r) cos(- 
B±(r) sin(- 




(5) 



and a uniform axial field 




(6) 



Here (r, <f>, z) are the cylindrical coordinates. The mag- 
nitude B±(r) is proportional to r near the axis r ~ 0; 
B±(r) ~ B'r, with B' being a constant. The system is 
assumed to be uniform along the z direction for calcula- 
tional simplicity. Our analysis should apply to a cigar- 
shaped system as well. 

Suppose we reverse the B z field slowly, while keeping 
B± fixed as shown in Fig. 1. It was shown in [5, 6] that 
a vortex with two units of circulation will be formed if we 
start with a vortex-free BEC. In these papers, an optical 
plug was introduced along the vortex axis to prevent the 
atoms from escaping from the trap when B = B± + 
B z vanishes at r = 0,< = T/2. Accordingly, the order 
parameter remains within the weak-field seeking state 
(WFSS) throughout the scenario. 

We suspect, however, that the BEC may be stable even 
without the optical plug since B vanishes only in r = 
at the time t = T/2. This will be justified by solving the 
Gross-Pitaevskii equation numerically below. In contrast 
with the previous work, we have to take the full degrees of 
freedom of the order parameter into account since the en- 
ergies of the three hyperfme states (the strong-field seek- 
ing state, the weak-field seeking state and the neutral 
state) are degenerate when B = 0. 
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FIG. 2: Polar angles (a,0) parametrizing B and I || F. 



B. Initial State 

We first solve the Gross-Pitaevskii equation in the sta- 
tionary state to find the initial order parameter configu- 
ration and the corresponding chemical potential. There 
is a sufficient gap between the weak-field seeking state 
and the other states at t = and the BEC may be as- 
sumed to be purely in the weak-field seeking state. Let 
us parametrize the external magnetic field B with use of 
the polar angles a and (3 as 



sin (3 cos a 

B(r,t) = B x {r) +B z (t) = \B\ | sin/3sina 

cos/? 

where a — — <p for the quadrupole field and 

\Bx(r)\- 



(3 = tan 1 



\B z {t)\ 



(7) 



(8) 



Since the hyperfine spin F is antiparallel with B, we 
must have 



— sin (3 cos a 
I = | — sin (3 sin a 
— cos/3 



(9) 



where I is a unit vector parallel to F, see Fig. 2. 

The most general form for the order parameter yielding 
the above I vector is 



, x , fo - 



cos (3 cos a + i sin a 
' I cos [3 sin a — i cos a 
— sin/3 



= /o««, (10) 



where {v a } represents the "phase", while fo is the am- 
plitude of the order parameter. In mathematical terms, 
{v a } defines the local 5*0(3) frame or the "triad" of the 
real orthonormal vectors {m, h, I}, where 



1 



in) a> I — rh x h. 



(11) 



It should be noted that the above order parameter takes 
the same form as that of the orbital part of the superfluid 
3 He-A. 

One may obtain more insight if the order parameter in 
Eq. ( |l0| ) is rewritten in the ^ basis as [13] 



9i 



foe 



-I'y 



1 ia 
2 e 



l (l - COS/3) 

^sin/3 
(1 + cos 13) 



(12) 



The angle (3 vanishes at r = (note that B± — at r = 
0), hence we find * +i = * = while oc e l ( Q -T) ^ 
there. For the order parameter to be smooth at r = 0, 
we have to choose 

7 = a = —(f). (13) 
Accordingly, the order parameter takes the form 



° V2 



in the ^ a basis and 



2 



cos f3 cos 4> — i sin </> 
- cos j3 sin (j) — i cos < 
— sin/? 



(1 - cos/3)e 2 ^ 
-V2sin/3e^ 
1 + cos (3 



(14) 



(15) 



in the basis. It follows from Eq. (|15J) that the state 
with I = —z {[3 = 0) has a vanishing winding number, 
while that with I = z ((3 = tt) has the winding number 
2. Accordingly, the former state may be continuously 
deformed to the latter state by changing (3 from to 
7r smoothly, resulting in the formation of a vortex with 
winding number 2 [14, 15]. It is shown in the Appendix 
that this phase may also be understood in terms of the 
Berry phase associated with the adiabatic change of the 
local magnetic field. 

The stationary Gross-Pitaevskii equation is given by 



2m 



ie abc ^ b uj Lc = j2^ a . 

(16) 



By substituting the WFSS order parameter in Eq. (|Tj) 
into Eq. (|l|), we obtain 



2m 



2 "?- A (7 - 8 cos /3 + cos 2/3) / 
2 4r z 

+lK B ( r )fo+9if6 i = A*/o- 



(17) 



Note that there appears an extra term f3' 2 that is missing 
if this BEC is described by a scalar (namely U(l)) order 
parameter. This term originates from v\V 2 v a and repre- 
sents the rotation of the local 50(3) frame. In ordinary 
physical settings, however, this term is smaller than the 
B-term by a factor of uj/ujl ~ 10~ 3 and hence its effect 
may be negligible in the following arguments. 
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The amplitude of the external magnetic field has an 
approximately harmonic potential profile near the origin, 



B{r,t) = yfB 2 ± (r) + B 2 (t) ~ y/B> 2 r 2 + B 2 (t) 
B' 2 , 



B,(t) 



2B z (t) 



(18) 



It turns out to be convenient to scale the energies by the 
energy-level spacing fiw at t = and the lengths by the 
harmonic-oscillator length oho, where 



= B' 



7/* 



mB z (0y 



oho 



(19) 
(20) 



If we take 23 Na and substitute B z (0) = 1 G and B' = 
300 G/cm, we find ajjo ~ 9.14 x 10 _1 fim and fiw ~ 
3.49 x 10 - 24 erg. For the same choice of the parameters, 
we have w L (r = 0,t = 0) ~ 4.40 x 10 6 rad/s - 1.330 x 
10 3 u. It is reasonable to assume r = 2-K/u)i,{r — 0,t = 
0) ~ 1.43 fj,s to be the measure of the adiabaticity. 

After these scalings, the dimensionless Gross- 
Pitaevskii equation takes a simpler form 



1. 



;V 2 / + 



B(r) 



P' 2 1 



! cos /? + cos 2/3) 



-5i Jo = fifo 



fo 
(21) 



where f = r/a H o,/o = /oa H 0'A = m/^, 5 = -f^B/hu 
and <h = gi/a^ huj ~ 0.0378. A similar scaling for 172 
yields 32 ~ —0.00132. Note that the singularity at f = 
vanishes if (3(f) approaches to fast enough as r — ► 0. 
The tildes will be dropped hereafter whenever it does 
not cause confusion. The eigenvalue /i may be obtained 
numerically. For 23 Na with B z and B' given above, we 
find [L — Hu>l(t — 0,t — 0) = 3.66, which amounts to 
1.28 x 10~ 23 erg in dimensional units. Figure 3 shows 
the corresponding condensate profile fo(r). 



C. Time Development 

The time development of the condensate wavefunction 
is obtained by solving the Gross-Pitaevskii Eq. (Q) which 
is written in dimensionless form as 



1 



V 2 * a +.g 1 |lf* Q 



+.g 2 (*) 2 ** +ie abc ^bUJLc 
The initial condition is 

*a = fo(r)v a , 



(22) 



(23) 



where fo(r) is obtained in the previous subsection and 
v a is defined in Eq. (not) . The condensate cannot remain 




FIG. 3: Condensate wave function fo{r). The condensate 
wave function and the radial coordinate are scaled by 
and oho, respectively, and hence dimensionless. 



within the weak-field seeking state during the formation 
of a vortex and we have to utilize the full degrees of free- 
dom of the order parameter ^> a . The Gross-Pitaevskii 
equation has been solved numerically to find the tempo- 
ral development of the order parameter while B z changes 
as 



B z (t) = B z (0) 



1 



2t 
f 



(0 < t < T), 



(24) 



with several choices for the reversing time T. Adiabatic- 
ity is not guaranteed at t ~ T/2 when the energy gaps 
among the weak- field seeking state (WFSS), the neutral 
state (NS) and the strong-field seeking state (SFSS) dis- 
appear. We expect, however, that this breakdown of adi- 
abaticity is not very significant to the condensate since 
it takes place only along the condensate axis for a short 
period of time. 

We project * a thus obtained to the local WFSS, NS 
and SFSS by defining the projection operators to the 
respective states as 



n Wa b=u Ufe ) II; 



Sab 



U a U, 



n 



Nab 



= S a b — Ilwab — lis 



ab 



where 



1 , „ ... 1 , A ... 

v a = —={m + in) a ,u a = — ={m - in) a 



V2' 



V2 K 



(25) 



(26) 



define the local WFSS and SFSS SO (3) frames, respec- 
tively. 

Figures 4, 6, and 8 depict the amplitudes of the 
wave functions at t = 0, T/2 ~ r, T/2 + r and T for 
T = 10r, 100r, and lOOOr. Figures 5, 7, and 9 show the 
particle numbers within the trap for the same choices of 
T . To simulate the trap loss, we introduced a function 



h(r) 



1 — tanh 



r - r 
A 



(27) 



with r = 30 and A = 2. The wave function $f a (r) is 
multiplied by h(r) after each step of the Crank-Nicholson 
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FIG. 4: Temporal development of the order parameters 
|*wfss|, |^ns| and |*sfss| for T = lOr, where r ~ 1.43ps. 
The graphs are plotted for t = (initial state), t — T/2 — r 
(slightly before B z — is crossed), t = T/2 + t (slightly after 
B z = is crossed) and t = T (B z is completely reversed). 
Furthermore, B z is fixed to B Z (T) = -B z (0) for t > T. The 
bottom row shows the order parameters at t = 1200T. The or- 
der parameters and the radial coordinate are plotted in units 
of a HO and sho, respectively. 



algorithm. The parameter ro roughly corresponds to the 
trap size, while A is taken large enough to prevent the 
wavefunction being reflected at r ~ ro. Thus particles 
that reach beyond r ~ ro disappear from the trap. The 
parameters ro and A are introduced for purely compu- 
tational purposes and should not be confused with any 
realistic experimental settings. As t ^> T, the conden- 
sates in SFSS and NS disappear from the trap and the 
particle number reaches its equilibrium value. The bot- 
tom row of Figs. 4, 6, and 8 shows the wavefunctions 
after the equilibrium is attained. 

Figure 4 shows the behavior of the condensate when 
T = lOr. In this case, B z (t) is reversed so fast that most 
of the particles remain in the WFSS at t — T/2 — t and 



1-0 h 
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FIG. 5: The ratio N(t)/N(0) as a function of t/T for the re- 
versing time T — lOr. Here N(t) is the number of atoms/unit 
length along the vortex axis at time t. The magnetic field B z 
decreases monotonically for < t/T < 1, but is kept fixed to 
-B z (0) for T < t, see Fig. 1. 
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FIG. 6: Same as in Fig. 4, but for T = lOOr, r ~ 1.43/Lts. The 
graphs are plotted for t = 0, t = T/2 -r,t = T/2 + r,t = T, 
and t = 50T. 
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FIG. 7: Same as in Fig. 5, but for T = lOOr. 



FIG. 9: Same as in Fig. 5, but for T = lOOOr. 
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they are suddenly converted into the SFSS at t = T/2+t. 
These particles in the SFSS and the NS are eventually 
lost from the trap as t — > oo, see Fig. 5. 

Figures 6 and 8 show that the behavior of the conden- 
sate is qualitatively similar for T = 100t and T = lOOOr. 
The condensate is transferred from the WFSS to the 
SFSS more efficiently at t = T/2 - r for T = 1000t 
which leads to more WFSS components at t = T/2 + r 
in this case. At a later time t > T, the condensate is 
found to oscillate in the trap with the frequency ~ lo. 

We have also analyzed the case T = 10000t and found 
no qualitative difference compared to the case T = lOOOr. 
There are slightly more of the NS component in the for- 
mer case at t = T/2 — r than the latter case, which leads 
to less particles in its equilibrium state at t T . 

It is certainly desirable to have more particles remain- 
ing in the trap when a vortex is created. Figure 10 shows 
the ratio of the final particle number to the initial particle 
number as a function of the reversing time T. Note that 
the final particle number is evaluated when the equilib- 
rium is reached. It is found that a considerable amount 
of the condensate (> 1/3) is left in the trap for a wide 
range of the reversing times 10 2 < T/t < 10 4 . 

In summary our numerical calculations indicate that 
a large fraction of the BEC remains in the trap when 
B z is reversed with proper choices of the reversing time 
T. Although there remain SFSS and NS components at 
t = T, they eventually disappear from the trap at t 3> 
T when the equilibrium is reached. The condensate is 
converted into a state with two units of circulation. Thus 
we conclude that a vortex can be created, even in the 
absence of an optical plug, by simply reversing the axial 
field B z . We expect that this makes the experimental 
realization of the present scenario much easier than that 
with an optical plug [5, 6]. 



FIG. 8: Same as in Fig. 4, but for T = lOOOr, r = 1.43/iS. 

The graphs are plotted for t = 0, t = T/2 - r, t = T/2 + r, t= IV. CONCLUSIONS AND DISCUSSIONS 

T, and t = AT. 

We have analyzed a simple scenario of vortex forma- 
tion in a spinor BEC. An axial magnetic field is slowly 
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FIG. 10: Ratio of the residual particle number to the initial 
particle number as a function of log 10 (T/t). The crosses are 
the results of our numerical calculation while the interpolating 
curve is introduced as a guide. 

reversed while the quadrupole field has been kept fixed 
in a Ioffe-Pritchard trap, which results in a formation 
of a vortex with the winding number 2. The spinor de- 
grees of freedom have been fully utilized which renders 
our scheme much simpler compared with other methods. 

Since the vortex thus created has a higher winding 
number, it is metastable and eventually breaks up to two 
singly quantized vortices. The lifetime of the metastable 
state is quite an interesting quantity to evaluate since its 
magnitude, compared with the field-reversal time T and 
the trap time, will be crucial to determine the ultimate 
fate of the vortex. 

A similar analysis with an F = 2 BEC is under 
progress and will be published elsewhere. 
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APPENDIX A: VORTEX FORMATION AND 
BERRY PHASE 

It is shown that the formation of a vortex in our sce- 
nario is understood from a slightly different viewpoint. 
Namely, we show that the phase appearing in the end of 




F z =-l (r=0) 
(a) (b) 

FIG. 11: Formation of a vortex in the present scenario may be 
understood in terms of the Berry phase associated with each 
spin. When \B Z \ S> \B± \ at t = 0, all the spins near the origin 
have F z ~ — 1 with the same phase. Due to the "twisting" of 
the spins during the evolution, they obtain different phases 
depending on the trajectories they follow, (b) shows four 
spins at t — T/2 where B z vanishes. The spin 1 follows the 
trajectory denoted by Ci in (a) starting from the south pole 
and ending at the north pole. Similarly the spins 2, 3 and 4 in 
(b) follows the trajectories C2,Cz and d in (a), respectively. 
The shaded area (7r/2 < <j> < tt, < /? < n) in (a) is the solid 
angle subtended by the paths Ci and C2. Since this area 
measures n, the spin 2 has a phase tt relative to that of the 
spin 1 at t = T. Similarly the spins 3 and 4 obtain phases 2n 
and 37r, respectively, relative to that of the spin 1 at t = T. 



the process may be identified with the Berry phase asso- 
ciated with the adiabatic change of the magnetic field. 

Let a point on the unit sphere in Fig. 11 (a) denote 
the hyperfme spin state with F = 1. All the spins near 
the origin r = have F z ~ — 1 at t = and hence corre- 
spond to the points near the south pole. Now the axial 
field B z monotonically decreases such that eventually all 
the spins near the origin have F z ~ +1 at t — T and 
hence are expressed by the points near the north pole. 
The path a spin follows on the unit sphere depends on 
the position of the spin relative to the origin. Figure 
11 (b) shows the configuration of four hyperfme spins F 
at t = T/2, when B z vanishes and F is determined by 
the quadrupole field. The path C\ in Fig. 11 (a) shows 
the trajectory the spin 1 in Fig. 11 (b) follows when t 
is varied from (/3 = 0) to T (/3 = 7t). Similarly the 
path C2 in Fig. 11 (a) is the trajectory of the spin 2 in 
Fig. 11 (b), and so on. When the spins 1 and 2 left the 
south pole at t = 0, they had the common phase factor, 
while at t = T they obtain the relative phase equal to 
the solid angle subtended by the trajectories C\ and Ci. 
The shaded area in Fig. 11 (a) shows this area which 
subtends the solid angle tt. Similarly the paths 1 and 3 
(1 and 4) subtend the solid angle 2-7T (37r), resulting in the 
relative phase 27r(37r) between the spins 1 and 3 (1 and 
4), respectively, when B z is completely reversed. Accord- 
ingly, as one completes a loop surrounding the origin, one 
measures the phase change of 47r observing that a vortex 
formation has taken place. 
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